A minimal two-band model for the superconducting Fe-pnictides 
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Following the discovery of the Fe-pnictide superconductors, LDA band structure calculations 
showed that the dominant contributions to the spectral weight near the Fermi energy came from 
the Fe 3d orbitals. The Fermi surface is characterized by two hole surfaces around the F point and 
two electron surfaces around the M point of the 2 Fe/cell Brillouin zone. Here, we describe a 2-band 
model that reproduces the topology of the LDA Fermi surface and exhibits both ferromagnetic and 
q — (ir, 0) spin density wave (SDW) fluctuations. We argue that this minimal model contains the 
essential low energy physics of these materials. 

PACS numbers: 71.10.Fd, 71.18.+y, 71.20.-b, 74.20.-z, 74.20.Mn, 74.25. Ha, 75.30.Fv 



Introduction - The recent discovery of superconductiv- 
ity in a family of Fe-based oxypnictides with large tran- 
sition temperatures [H, 0, d, 0, [B|, Q has led to tremen- 
dous activity aimed at identifying the mechanism for su- 
perconductivity in these materials. Preliminary exper- 
imental results including specific heat [a, point- contact 
spectroscopy [1] and high-field resistivity [jllld] measure- 
ments suggest the existence of unconventional supercon- 
ductivity in these materials. Furthermore, transport 11 
and neutron scattering [lij] measurements have shown the 
evidence of magnetic order below T — 150K. An experi- 
mental determination of the orbital and spin state of the 
Cooper pairs, however, has not yet been made. 

The high transition temperatures and the electronic 
structure of the Fe-pnictide superconductors suggest that 
the pairing interaction is of electronic orig i n. II 311 First- 
principle band structure calculations LJ, lla . I la . Il7j ] have 
shown that superconductivity in these materials is asso- 
ciated with the Fe-pnictide layer, and that the density of 
states (DOS) near the Fermi level gets its maximum con- 
tribution from the Fe-3d orbitals. The consensus based 
on these calculations is that the Fermi surface consists of 
two hole pockets and two electron pockets. Calculations 
from Ref. fl5i | also show van Hove singularities which 
might be responsible for enhanced ferromagnetic fluctu- 
ations. The bare magnetic spin susceptibility determined 
from these bands exhibits both ferromagnetic q ~ and 
finite q SDW peaks. 

Several tight binding models for the band structure 
have been proposed. Cao et al. [HI used 16 localized 
Wannier functions to construct a tight-binding effective 
Hamiltonian. Kuroki et al. [3] have used a 5 orbital 
tight binding model to fit the band structure near the 
Fermi ener gy. Others have introduced generic 2-band 



models [2fJ, |2l|, |22j. However, the relationship of these 



latter models to the multiple Fermi surface electron and 
hole pockets found in LDA calculations is unclear. Since 
it appears likely that these multiple Fermi surfaces play 
an essential role in determining the momentum depen- 



dence of the spin and orbital fluctuations which would 
mediate an electronic pairing mechanism, we would like 
to construct a minimal model that exhibits a Fermi sur- 
face similar to that obtained from band structure calcu- 
lations. 

This model has two orbitals per site on a two dimen- 
sional square lattice. By adjusting the one-electron hop- 
ping parameters and the chemical potential one can ob- 
tain a Fermi surface which has the same topology as 
found from the band structure calculations. The non- 
interacting spin susceptibility also exhibits both ferro- 
magnetic and finite q SDW peaks. With the addition of 
an onsite intra-orbital and inter-orbital Coulomb inter- 
actions, and an intra-orbital Hunds rule coupling, this 
model represents what we believe is a minimal model 
for describing the low energy physics of these materi- 
als. In addition, the relative simplicity of this model 
should be useful in the phenomenological analysis of ex- 
periments related to the gap symmetry [23| and in nu- 
merical density-matrix renormalization group (DMRG) 
and dynamic cluster studies. 

Model Hamiltonian - The structure of the FeAs layer 
of LaOFeAs viewed along the c-axis is illustrated in 
Fig. la. The Fe ions form a square lattice which is inter- 
laced with a second square lattice of As ions. These As 
ions sit in the center of each square plaquette of the Fe 
lattice and are displaced above and below the plane of 
the Fe ions as indicated in the figure. This leads to two 
distinct Fe sites and a crystallographic unit cell which 
contains two Fe and two As ions. As shown by various 
band structure calculations, the main contribution to the 
density of states within several eV of the Fermi surface 
comes from the Fe 3d states which disperse only weakly 
in the z-direction. The 3d Fe orbitals hybridize among 
themselves and through the As p orbitals leading to a 
complex of bands. However, as noted in Ref. [161 ]. the 
bandstructure near the Fermi level is relatively simple in 
the unfolded lFc/ccll BZ where it primarily involves three 
Fe orbitals d xz , d yz and d xy (or d x 2_ y 2). Based upon this 
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FIG. 1: (a) The Fe ions form a square lattice and the crys- 
tallographic unit cell contains two Fe and two As ions. The 
As ions are located either directly above (solid circles) or be- 
low (dashed circles) the faces of the Fe square array. (b)A 
schematic showing the hopping parameters of the two-orbital 



d yz model on a square lattice. Here ti is a near neigh- 



bor hopping between <r-orbitals and £2 is a near neighbor hop- 
ping between 7r-orbitals. We also include a second-neighbor 
hopping £4 between different orbitals and a second-neighbor 
hopping £3 between similar orbits. 



observation and by making the further approximation 
that the role of the d xy (d x 2_ y 2) orbit can be replaced by 
a next near neighbor hybridization between d xz , d zy or- 
bitals, we consider a two-dimensional square lattice with 
two degenerate "d^z, dyz" orbitals per site. While one 
may well need a third orbit to control the relative sizes 
and eccentricities of the electron and hole pockets, we 
find that a two-orbital model can lead to a Fermi sur- 
face which resembles that obtained in the bandstructure 
calculations. 

The tight-binding parameters of the 2-orbital model 
that we will study are illustrated in Fig. [TJ It is conve- 
nient to introduce a two-component field 



d xs (k) 
d ys (k) 



(1) 



Here d xs {k) (d ys (k)) destroys a d xz (d yz ) electron with 
spin s and wave vector k. Then the tight binding part of 
the Hamiltonian can be written as 



H = ^2 ^ts i( £ +( k ) - l«) 1 + £-(k)r 3 + e xv {k)Tx] tp ks 



k S 



(2) 



with Ti the Pauli matrices and 
e x (k) ± s y (k) 



s±(k) 
£x{k) 

Ey{k) 

£xy(k) 



-2t\ cos k x — 2ti cos k y — cos k x cos k y , 
-2^2 cos k x — 2t\ cos k y — cos k x cos k y , 
-4t4 sin k x sin k v . 
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FIG. 2: (a) The band structure of the two-band model with 
£1 = — 1, <2 = 1-3, £3 — £4 = —0.85 and fi = 1.45, plotted along 
the path (0, 0) -» (it, 0) -> (tt, tt) -> (0, 0) as shown in Fig. |3] 
(a) by the black dashed lines, (b) The band structure folded 
to the small BZ, with the F, X, M defined in the small BZ as 
shown in Fig. 02(b). (c) The two-d band structure for k x ,k y £ 
[0, tt]. A saddle point exists for each band, (d) The density of 
states of the two band model, with two Van Hove singularities. 
The dashed line shows the fermi level corresponding to our 
choice of fi = 1.45. 



The one electron Matsubara Green's function is given 
by 

G.(fe,i"n) = {lUJn , - e+ik)) |,T " (fc) " 3 = ?A (fc)fl (3) 



(iuj n - E+(k)) [iun - E-(k)) 



with 



E ± {k) = e+(fc) ± Jel(k) + eUk) - fi 



(4) 



In Figure we show the band structure of the model 
for a specific choice of hopping parameters t\ = —1, £2 = 
1.3, £3 = £4 = —0.85, in units of |£i|. The folded energy 
spectrum in Fig. [|] (b) shows the band structure in the 2 
Fe/cell zone. Due to the saddle points in the energy spec- 
trum (as shown in Fig. [U (c)), there are two Van Hove 
singularities in the density of states, which also qualita- 
tively agrees with the LDA results. 15j In Figure [3] we 
show the Fermi surface for the same set of parameters. 
On the large BZ (Fig. [3k,) associated with our model 
which has 1 Fe/unit cell, there are two hole Fermi pock- 
ets labeled a.\ and a 2 defined by E-(k) — 0, and two 
electron Fermi pockets (3\ and (3% defined by E + (k) = 0. 
To compare with band structure calculations, one must 
fold the large BZ into a smaller one which is dual to the 
crystallographic unit cell containing two Fe atoms. The 
dashed square in Fig. [3] (a) marks this smaller zone and 
in Fig. [3] (b) we show what happens as the and ftx,2 
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be obtained as 

T 

Xs(q,*fi) = X/ Tr[G(fe + g,io;„+iO)G(fe,iw n )] 



(k + q, i/|k, ^ 



/\|2 



2N ifl + £Wq - £v,k 

k,ZAi/' 

(^(-Br/.k+q) - ra-F k)) 



(6) 



Here E^k, " = 1) is the eigenvalue of the upper 

(lower) band given by Eq. and |k, v) the correspond- 
ing eigenvector. tif(E) = l/(e l3E + 1) is the fermi distri- 
bution function. 
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FIG. 3: (a) The Fermi surface of the 2-orbital model on the 
large lFe/cell BZ. Here, 01,2 surfaces are hole Fermi pock- 
ets given by E-(kf) = and /3i,2 are electron Fermi pockets 
by E+(kf) = 0. The dashed square indicates the BZ of the 
2Fe/cell. (b) The Fermi surface folded down into the 2 Fe/cell 
BZ consists of two a surfaces around F and two elliptically de- 
formed j3 surfaces around the M point. Here the parameters 
are the same as in Fig. [5] 



bands of Fig. [3] (a) are folded back into the 2Fe/cell BZ. 
One sees that this gives Fermi surfaces with the same 
topology that is obtained from LDA band structure cal- 
culations [24j. 

One-loop spin susceptibility.- Now we study the one- 
loop spin-susceptibility for the tight-binding model @. 
Due to the existence of two degenerate orbitals in our 
model, the spin susceptibility also has orbital indices, 
and is defined by 



FIG. 4: (a) The uj = bare spin susceptibility xs(q) versus 
q for the same tight-binding parameters as used in Fig. [3] 
(b) The bare spin susceptibility xs(q) (red solid line) and the 
RPA spin susceptibility Xs PA (q) for U = V = 6, J = (blue 
dashed line) along the (0, 0) — ► (tt, 0) — > (tt, n) —> (0, 0) path 
in BZ, as shown by the dashed line in subfigure (a). 

Fig. |4] shows a plot of the static spin susceptibility 
Xs(<Zi0) versus q, where one can see the structure as- 
sociated with the various nesting points and density of 
states features. For our choice of parameters, the largest 
value of xo(q) occurs around q — (tt, 0) and (0, tt), which 
suggests a transition to an antiferromagnetic (AFM) or- 
dered phase at some critical interaction strength. This 
is also in a gre ement with the result of band structure 
calculations 



1( 



Xs*(q, «n) 



p 



dre inT (r r S s (-q,r).S t (q,0)} (5) 



n 



251 ] . A recent neutron scattering exper- 
iment has confirmed that such a peak develops below 
T - 150K 0. Such a peak in the spin susceptibility 
comes from the nesting between the electron and hole 
Fermi pockets, which can be seen from the chemical po- 
tential dependence of the spin susceptibility. As shown in 
Fig. O the spin susceptibility at q = (0, 0) jumps discon- 
tinuously at /j, ~ 1.2, which follows the behavior of the 
density of states shown in Fig. [2] (d), and corresponds to 
the onset of electron Fermi pockets. At the same time, 
the q = (0, tt) spin susceptibility is also enhanced sig- 
nificantly due to the nesting effect. When the chemical 
potential is increased further, the fermi level gets closer 
to the Van Hove singularity, and the hot point of the spin 
susceptibility is shifted gradually to the neighborhood of 
(0, 0) and (tt, tt), as shown in Fig. [5] (b). 
s >* The RPA spin susceptibility. — Now we consider the 

effect of electron-electron interaction in this model. For 
The one loop contribution to the spin susceptibility can the two d orbitals we considered, the generic form of the 



here s,t = 1,2 label the orbital indices, and S s (q) = 
2 Ylk "^iL(k + 0)^/3^80 (k) is the spin operator for the 



orbital labeled by s. 
given by 



The physical spin susceptibility is 



xs(q,«n) =^x st (q,«n). 



4 




FIG. 5: (a) The chemical potential dependence of Rex(q, w = 
0) with q = (0,0), (0, tt) and (71", tt). The vertical dashed line 
shows the chemical potential /1 = 1.45 which we are work- 
ing on. (b) The three-d color plot of spin susceptibility as a 
function of k x , k y and chemical potential fi. The hot points in 
BZ shifts from (n,0) to (0,0) and (n, n) with electron doping 
as n increases. The deepest red region shows the Van Hove 
singularities. 



on-site interaction can be written as 



Hi, 



Vnurii2 - JSu • S* 



with U and V the intraband and interband Coulomb re- 
pulsion, and J the Hunds rule coupling. For an isolated 
Fe atom, the intraband U and interband V are similar 
in magnitude, and J is an order of magnitude smaller 
[26j - Thus we expect the U and V to be the dominant 
terms in the interaction. We suggest that H = H + H mt 
represents a minimal model for the Fe-pnictides super- 
conductors. 

Next we will study the effect of such interactions on 
the spin fluctuations within RPA. Due to the two band 
nature of the model we considered, the RPA correction 
should be calculated for the generic spin susceptibility 
Xst(q, iSl) defined in Eq. (O, which is determined by the 
following matrix equation: 



x KFA (q, in) = Xo (q, in) (i - r Xo (q, in)) - 



(8) 



Here xo is the 2x2 matrix formed by the intra-orbital 
and inter-orbital spin susceptibility defined in Eq. 
and r is the interaction vertex defined by 



U J/2 
J/2 U 



(9) 



We note that the interband interaction V does not con- 
tribute to the RPA response when only the spin fluctua- 
tions are considered. 

In the following we set J = 0, which makes the inter- 
action vertex T in Eq. ([9]) proportional to the identity. 
For the tight-binding model parameters used in Fig. [3] 
and U = 6, we obtain the RPA spin susceptibility shown 
in Fig. [J] (b) by the dashed line. As expected, the spin 
susceptibility is enhanced around the hot points (tt, 0) 
and (0, tt). We have also carried out the RPA calculation 



for a finite Hunds rule coupling J > 0, and find that the 
spin fluctuations are enhanced by increasing J, but the 
structure of x(q) remains qualitatively the same. 

In conclusion, we have described a minimal model for 
the Fe-pnictides which we believe contains the essential 
low energy physics of these materials. This model con- 
sists of a two dimensional square lattice of sites with 
each site having two degenerate orbitals. By fitting the 
tight binding parameters, one can obtain a band struc- 
ture which, after folded to the 2Fe/cell BZ, exhibits two 
hole pockets around the Y point and two electron pock- 
ets around the M point. The electron-electron interac- 
tions are taken to be onsite intra-orbital and inter-orbital 
Coulomb interactions U and V and an onsite Hund's rule 
coupling J. The structure of the bare spin susceptibility 
is peaked around (tt, 0) for the parameters we chose to 
fit the fermi surface. Such AFM spin fluctuations also 
leads to the possibility of non-conventional superconduc- 
tivity, which we will discuss in a separate work. Different 
types of spin or orbial orders and superconductivity can 
possibly occur for different fillings. Therefore, we con- 
clude that this model contains a rich variety of magnetic, 
orbital and pairing correlations. 
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